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基于 小 波 包 阐 值 法 的 脉冲 星 信号 去 噪 方法 研究 
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摘要 : 小 波 变换 去 噪 只 对 信号 的 低频 部 分 进行 分 解 和 处 理 ， 忽 略 了 信和 号 的 高 频 部 分 的 分 解 与 处 理 ， 无 法 对 
信号 的 整个 频段 信息 的 有 效 提取 ， 进 而 影响 脉冲 星 信号 去 噪 的 效果 。 鉴 于 此 ， 本 文 提出 利用 小 波 包 闪 值 法 
对 脉冲 信号 进行 分 析 ， 针 对 小 波 包 分 解 后 的 系数 选用 不 同 的 阀 值 函 数 处 理 。 实 验 结 果 表 明 ， 使 用 该 方法 进 
行 脉冲 星 信号 的 去 品 效 果 得 到 了 提高 ， 其 中 脉冲 星 信号 的 信 品 比 、 峰 值 信 噪 与 平滑 度 指标 都 有 相应 的 得 到 
了 改善 ， 为 脉冲 信号 的 去 噪 提供 了 一 种 新 的 思路 。 
关键 词 : 小 波 包 分 解 ; 国 值 法 ; 脉冲 星 信 号 ; XU 
中 图 法 分 类 号 P162.5 文献 标志 码 A 

2016 年 脉冲 星 观 测 研 究 利 器 FAST 在 贵州 平 塘 落成 ， 自 此 我 国 开始 了 依靠 国内 设备 发 现 
脉冲 星 的 零 突破 。 脉冲 星 是 宇宙 中 的 一 类 奇妙 天 体 ， 是 验证 强 引 力 场 、 强 磁场 和 高 密度 等 极 
端 物理 环境 下 物理 规律 的 “天 然 实 验 室 ”" 。 作 为 脉冲 星 研 究 的 利器 ，FAST 采用 19 波束 多 
波束 接收 机 系统 每 天 8 小 时 对 天 区 进行 大 规模 巡天 搜索 , 搜寻 并 接收 脉冲 星 发 出 的 电磁 波 脉 
冲 信号 , 依据 电磁 脉冲 信号 确定 天 区 中 是 否 存 在 脉冲 星 。 众 所 周知 ， 处 在 遥远 星空 中 的 脉冲 
星 距 地 球 约 3000”55000 光 年 , 脉冲 星 发 出 的 脉冲 信号 到 达 地 球 且 成 功 被 专用 仪器 接收 的 信 
号 能 量 十 分 微弱 ， 整 个 脉冲 信号 淹没 在 极 强 的 噪声 中 *”。 因 此 ， 脉 冲 星 信 号 的 去 噪 在 整个 
脉冲 星 信 号 的 检测 和 发 现 中 显得 至 关 重 要 。 
小 波 变换 作为 一 种 多 分 辩 分 析 方 法 能 够 同时 在 时 域 和 频 域 上 对 信和 号 进行 分 析 , 由 于 其 具 
有 时 域 局 部 化 和 多 分 辨 特 性 而 被 业内 人 士 称 之 为 “数学 显微镜 ”, 非常 适合 处 理 非 平稳 信和 号 
“。 脉 冲 星 信号 是 一 种 典型 的 非 平稳 信号 ， 朱 晓 明 “等 人 提出 小 波 分 析 运 用 到 脉冲 星 信号 去 
品 。 他们 认为 小 波 分 析 可 以 更 好 的 对 非 平稳 信号 进行 消 噪 处 理 , 因而 小 波 变换 蔡 代 了 传统 的 
傅 里 时 变换 , 在 脉冲 星 信号 的 处 理 上 大 显 身手 。 此 后 研究 者 们 对 小 波 变换 去 除 脉冲 星 信 和 号 噪 
声 进行 了 大 量 的 研究 ， 阁 迪 ™ 等 人 探索 了 极 大 极 小 原理 、Stein 无 偏 似 然 原理 的 脉冲 星 降 品 
方法 ; 柯 巾 政 " 等 人 尝试 了 改进 闵 值 函数 以 提高 脉冲 星 信 号 去 噪 效 果 ， 提 出 了 软 硬 折 中 准 值 
函数 、 指 数 型 闪 值 函数 处 理 脉冲 星 小 波 分 解 系数 ， 获 得 了 较 好 的 效果 。 

小 波 变换 既 能 够 较 好 的 去 除 信和 号 的 噪声 ,又 能 较 好 的 保持 原 有 信和 号 的 数据 , 是 处 理 信和 号 
的 相对 理想 的 工具 。 但 小 波 变换 主要 是 针对 低频 段 信息 的 提取 和 分 析 , 对 于 信号 的 高 频段 信 
息 则 经 常 忽略 掉 , 这 对 于 整个 信号 分 析 的 精确 度 产生 较 大 的 影响 。 小 波 变换 去 噪 是 在 小 波 函 
数 中 选择 一 种 合适 的 信号 处 理 基 沙 数 , 利用 小 波 基 函数 进行 小 波 分 解 , 根据 不 同 阔 值 函数 的 
处 理 结果 进行 对 信号 的 重 构 ， 从 而 达到 小 波 变换 去 品 。 小波 变 换 去 噪 在 小 波 基 函数 的 选择 和 
闵 值 函数 选择 的 不 同 而 出 现 不 同 的 去 噪 效果 ， 使 得 小 波 去 噪 效果 存在 了 一 定 的 随机 性 ”。 

随 着 小 波 理 论 的 进一步 发 展 , 小 波 包 分 析 逐 步 进 入 人 研究 者 们 的 视野 。 小波 包 分 析 是 小 波 
分 解 的 一 般 化 , 具有 更 好 的 信和 号 去 噪 能 力 ， 既 能 够 像 小 波 分 析 一 样 对 信和 号 的 低频 部 分 进行 分 
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解 ， 也 能 够 对 信号 的 高 频 部 分 进行 分 解 ， 能 更 好 地 提取 各 频段 的 有 用 信息 。 小 波 包 理 去 噪 方 
法 在 图 像 、 地 震 信号 、 心 电信 号 、 海 波 信 号 等 领域 不 断 应 用 ， 刘 秀平 等 "开始 尝试 将 小 波 包 
滤波 方法 与 X 射 线 脉冲 星 脉冲 噪声 滤 除 进行 结合 ， 但 并 没有 对 小 波 包 去 噪 进行 系统 的 分 析 ， 
鉴于 此 , 本文 提 出 将 小 波 包 阔 值 法 应 用 到 脉冲 星 信号 的 去 噪 分析 中 ,系统 分 析 各 种 阔 值 函数 
获得 的 去 噪 效果 ,实验 显示 小 波 包 冰 值 法 与 传统 的 小 波 变换 法 相 比 ,脉冲 星 信号 的 去 噪 效 果 
得 到 了 明显 的 提高 ， 为 脉冲 星 信 号 的 研究 提供 一 些 思 路 。 


1 小 波 变 换 与 小 波 包 变换 


1.1 小 波 变换 方法 
观测 信号 可 以 看 作为 f(t)， 满 足 : 
f(t)=s(t)+ n(t) (1-1) 
其 中 f(t ) 为 接收 到 的 原始 信号 即 含 曲 信 号 ，n(t ) 是 方差 为 o? 的 噪声 ， 服 从 n(0，oa?*) 分 布 ， 
s(t) 为 需要 提取 的 去 噪 后 的 信号 ,一 维 含 品 信 号 f(t) 按 照 小 波 级 数 展开 可 以 得 到 1n 个 分 辩 率 
小 波 的 线性 组 合 ， 有 具体 的 小 波 变 换 为 : 


Wr Win) = 22 So Fo(2in-x) dy (1-2) 


© 其 中 Wi (wjx) 为 2(R) 中 取 到 的 小 波 基 册 jx 的 离散 小 波 变换 ， 小 波 变换 系数 由 上 式 确定 。 
= 1.2 小 波 包 变换 方法 
e 小 波 包 变换 是 小 波 变 换 的 拓展 , 具备 小 波 变换 处 理 观测 信号 低频 部 分 的 能 力 , 同时 能 够 
QN 较 好 的 分 解 观测 信号 的 高 频 部 分 ， 收 集 观测 信号 各 个 频段 的 有 用 信息 ， 进 而 提高 去 噪 精度 。 
o 具体 的 小 波 包 变换 为 : 
N P ara 
W(x) = V2 Ynez 8, 92x — n) 
其 中 op(CO 和 uCO 分 别 为 信号 在 LzCR) ELD P RE BAT) ERE 
数 为 各 和 g ， 他 们 分 别 满 足 的 条 件 为 ; 


(1-3) 


IL 


中 滤波 器 系 


Ynez hn (2x - n) = Ok (1-4) 
Enez hn = V2 (1-5) 


— BaT (—1)°hn-1 (1-6) 
2 小 波 包 多 阅 值 去 噪 ; 
2.1 小 波 包 分 解 与 重 构 


小 波 包 理论 中 的 小 波 包 变换 具有 分 解 与 重组 的 性 质 中 ， 具 体 小 波 包 分 解 算法 的 数学 表 
达 式 为 : 


WA (k) = Xiez hi- 2x wi (D (2-1) 

wrt (k) = Diez & 2, W O (2-2) 
Few) (表示 在 分 解 尺度 为 上 的 第 mn 条 分 又 树 上 的 第 上 个 小 波 包 分 解 系数 ， 而 ww 向 Ck) 
和 w 台 (各 表示 w? (Co) 的 两 个 分 又 ， 具 体 相应 的 小 波 包 重 构 算法 的 数学 表达 式 为 : 

wj (k) = Dez hy-ai we D + Xiez 8k-21 wes) (2-3) 
其 中 w 为 信号 进 过 小 波 包 分 解 后 得 到 的 分 解 系数 ，h、g 为 小 波 包 分 解 的 滤波 器 系数 ; 
j、7n 为 小 波 包 分 解 中 各 分 解 节点 的 编号 ;1、Kk 则 为 信号 进行 小 波 包 分 解 所 在 的 分 解 层 数 。 
2.2 最 优 小 波 包 基 的 确定 
小 波 包 分 析 处 理 观 测 信号 时 , 对 信和 号 进行 小 波 包 分 解 。 小 波 包 分 解 与 小 波 分 解 有 所 不 同 ， 

小 波 包 分 解 是 按照 二 又 树 的 形式 进行 分 解 , 信号 的 低频 和 高 频 部 分 都 进行 相应 的 完全 二 又 树 


的 形式 进行 分 解 。 信 和 号 分 解 的 层 数 会 直接 影响 信号 的 去 品 效 果 , 处 理 的 信号 在 进行 逐 层 


过 程 中 必须 要 确 


定 


个 最 优 分 解 层 数 。 分 解 层 次 过 低 会 引起 频段 中 混 有 有 用 信号 ,分 解 


分 解 


= 
2 


^ 


过 高 时 亦 会 使 得 结果 有 误差 5 。 在 小 波 分 析 去 噪 过 程 中 分 解 层 数 过 多 会 造成 信号 信息 量 的 


丢失 ， 从 而 导致 信 品 比 下 降 ， 影 响 最 终 的 计算 效果 。 在 小 波 分 析 的 实际 应 用 中 ， 分解 尺度 一 
般 的 取 3~502， 在 本 文 的 研究 中 取 作 5 进行 对 比分 析 。 


小 波 包 变换 9 
价 函 数 的 评价 标 ; 
给 定 区 间 的 紧 文 


真正 的 目 


HUY SymN. CoifN, DbN. Bior 系列 小 波 以 及 Haar 小 波 进行 了 比较 分 析 。 
析 ， 本 节 中 统一 采用 固定 阐 值 规则 ， 利 用 不 同 的 小 波 基 进行 5 层 分 解 ， 分 忆 


2、 表 3、 表 4 所 示 。 其 中 表示 小 波 阶 数 、SymN、CoifN、DbN、BiorN 分 别 表 示 Symlets、 


FP 存在 着 多 种 小 波 构 成 的 基底 库 , 因此 需要 根据 所 处 理 信号 的 特 怕 
住 进行 选择 。 选 择 小 波 函 数 时 需要 充分 考虑 脉冲 信号 的 特点 ,选择 能 够 满足 
性 和 是 够 的 消失 和 矩 ,以 便 能 够 更 加 有 效 二 
标 脉 冲 信 号 。 为 了 从 符合 脉冲 星 信号 特点 的 小 波 基 函数 中 硬 


Coiflets、Daubechies、Biorthogonal 系列 小 波 。 
表 2.1 SymN. CoifN, DON 系列 小 波 基 的 脉冲 星 信 号 去 噪 后 信 品 比 
Table 2.1 The SNR of Pulsar signal denoised by SymN, CoifN and DbN series wavelet bases 


Sym N Coif N DbN 
阶 数 N 

SNR SNR SNR 
1 04, 495421 93. 487519 
2 104. 028278 22. 461700 04. 028278 
3 113. 023734 137. 609002 113. 023734 
4 117. 713573 50. 993809 121. 279292 
5 128. 273024 163. 244416 128. 628087 
6 135. 525030 74. 287862 35. 497251 
7 142. 075414 176. 879722 142. 056956 
8 148. 166205 172. 385781 48. 189196 
9 153. 963343 171. 814606 154. 003296 
0 159. 679101 72. 627403 59. 637559 
11 161. 617158 175. 312769 165. 295176 
2 170. 280179 76. 374491 70. 247172 
13 173. 439558 180. 927850 172. 002886 
4 178. 031220 77. 493286 77. 478973 
15 176. 211857 176. 077450 182. 380212 


K 2.2 SymN. CoifN. DbN 系列 小 波 基 的 脉 六 
Table2.2 The SNR of Pulsar signal denoised by SymN, CoifN and DbN series wavelet bases 


Paes vere ia 


E， 结 合 代 


也 去 除 脉冲 信号 掺 杂 的 噪声 , 提取 出 
DE ABE PBL, 本 节选 
为 了 便于 比较 分 


FARUK 1、 表 


Sym N Coif N DbN 
Brae N 

RMSE RMSE RMSE 
1 1.1590260180636194e+38 1.159937964771207e+38 
2 1.159105 137508623e+38 1.1591729484203662e+38 1.159105 1375086237e+38 
3 1.15902630444004e+38 1.159201733565514e438 1.1590263044471583e+38 
4 1.1591290190598036e+38 1.1591803200172715e+38 1.1593325604448439e+38 
5 1.1589381807874201e+38 1.1591317586993424e438 1.159134635851061e+38 
6 1.1590973 14892064 1e+38 1.1590666107944755e+38 1.1590291941601505e+38 
yi 1.159204508875512e438 1.1590039338852231e438 1.1591823177722796e+38 
8 1.159088 1427106949e+38 1.1589613392166439e+38 1.1591410945254123e+38 


1.158947 


1.158964702074117e+38 


1.159004. 
1.159049 
1.159079. 
1.159075 


1.159035404037832e+38 


9802663968e+38 1.159008 1322083972e+38 


1.1590547974489843e+38 


2562137324e+38 1.1591076638631676e+38 
5391531836e+38 1.1590266700820342e+38 
8965814301e+38 1.1590067783780341e+38 
9265506637e+38 1.15905305165149e+38 


1.1590402208654303e+38 


9 1.1589676535028293e+38 
10 1.1590847694307965e+38 
11 1.1589516584680341e438 
12 1.1591001911149105e+38 
3 1.158990764381691e438 
14 1.1590776721200966e+38 
5 1.1590080833 152019e+38 
K 2.3 Bior N 系列 小 波 基 的 脉 六 


是 信号 去 噪 后 信 品 比 和 均 方 根 误差 


Table 2.3The SNR and RMSE of Pulsar signal denoised by SymN, CoifN and DDN series wavelet bases 


Bior N 
阶 数 入 
SNR RMSE 
1.1 93.487519 1.159937964771207e+38 
1.3 92.882619 1.161230305420237e+38 
15 92.564777 1.1621428245 153362e+38 
2.2 113.496536 1.1620257884498562e+38 
2.4 114.275923 1.1593361263215988e+38 
2.6 114.380095 1.1591729507026387e+38 
2.8 114.370066 1.15925 14421028786e+38 
3.1 121.556228 1.4888560080440952e+38 
3.3 31.610178 1.1712494267787645e+38 
3.5 133.428810 1.1611795518841145e+38 
347 33.791348 1.1596986873497866e+38 
3.9 133.895611 1.1592749724774778e+38 
44 28.219091 1.159259591601108e+38 
35.5 139.266704 1.1592072048905078e+38 
6.8 143.374262 1.15912528098707e+38 
X24 Her 系列 小 波 基 的 脉冲 星 信号 去 噪 后 信 品 比 和 均 方 根 误差 


Table2.4 The SNR and RMSE of Pulsar signal denoised by Haar wavelet bases 


Haar 


阶 数 N 


SNR 


RMSE 


93.487519 


1.159937964771207e+38 


从 以 上 四 张 表 中 统计 的 均 方 误 根 误差 (RMSE) 和 信 品 比 (SNR) 可 以 看 出 利用 Sym 
CoifN、DbN 、Bior N. Haar 5 种 小 波 进 行 脉冲 星 信 号 去 噪 的 整体 上 较 好 ， 其 中 SymN 系列 


小 波 的 效果 最 好 ，CoifN 小 波 次 2 
同 ， 单 个 小 波 基 而 言 ，Sym5 去 噪 效果 最 好 。 
Symmlet5 作为 小 波 基 函 数 。 因 此 对 脉 ; 


! 信 号 进 


小 波 包 基 ， 也 就 是 小 波 包 分 解 的 最 优化 树 。 


2.3 多 阅 值 函 数 的 选择 


利用 小 波 包 变换 处 理 信号 的 过 程 


经 常会 出 现 “ 笨 


，BiorN 小 波 最 差 。 此 外 小 波 基 阶 数 不 同 其 去 噪 效果 也 不 
因此 本 文中 选取 了 近似 对 称 的 紧 支 集 正 交 小 波 


ll 


行 小 波 包 分 析 ， 分 解 层 数 为 5 次 ， 得 到 最 优 


器 时 会 进行 的 一 种 频率 交大 现 象 , 从 而 导致 分 解 的 小 波 包 系数 
号 的 低频 部 分 从 小 到 大 排列 ， 信 号 的 高 频 部 分 由 
理 之 前 需要 先 对 小 波 包 变换 后 的 小 波 包 系数 进 4 


的 信号 处 理 阔 值 。 


T 


j 转 ”现象 , 也 就 是 在 信号 经 过 高 通 滤波 
8 现 频率 大 小 顺序 错位 。 即 信 
分 则 是 从 大 到 小 排列 。 因 


此 在 进行 小 波 包 阔 值 处 


了 排序 , 再 根据 排序 好 的 小 波 包 系数 选取 对 应 


小 波 包 分 析 中 沿用 的 闵 值 选取 准则 有 四 种 ， 分 别 为 sqtwolog《〈 固 定形 式 ) 阔 值 准则 、 
rigrsure《 自 适应 〉 阔 值 准则 、heursure( 启 发 式 ) 阐 值 准则 、minimaxi〔 极 大 极 小 值 》 阐 值 
准则 ， 本 文采 用 sqtwolog〔 国 定形 式 〉 阔 值 准则 对 信号 进行 处 理 。 对 信号 分 解 后 每 层 的 系数 
选 定 阔 值 以 后 ， 需 要 确定 对 信和 号 处 理 的 阔 值 函数 。 和 党 用 的 阔 值 函数 有 软 阔 值 函数 、 硬 阔 值 函 
数 。 其 中 硬 阔 值 函数 采用 的 是 将 所 有 绝对 值 不 大 于 闪 值 X 的 所 有 元 素 用 零 取代 ， 这 样 容易 造 
成 处 理 后 的 信号 不 连续 。 硬 阔 值 函数 表达 式 为 


UN ua Iw?| z A 


y ( _ 
0, Iw? | <A 


wi = (2-4) 


A BREL B CU) v JL T 8 DR EE ANAL, EE ORE i fia “STAT AE, BCE EAE HE. PA 
而 在 面临 小 波 系 数 较 大 时 , ZEEE PI T EJ A 3 Fr B rj PUE D ZA HE TT SEMIS 
fam EMI. XI BARCA IJ 
= sign(w;') Kw?| 20 , |w?| 2A 
w” = 

7 0, Iw? | <i 


(2-5) 


AE JC. WEBEL PS CAE a SIN AREER ASAE, HER EB TE, 研究 人 员 提 出 了 
软 硬 阔 值 折 中 法 ， 有 具体 的 表达 式 为 
"S sign(w;') CK]w?| 20 , |w*| z X 
" =lo, Iwr| «A 


(2-6) 


式 子 中 的 a 取 值 范围 为 0 < a < 1. MAEA 0 ATEREA, MEA 1 时 相当 
于 软 闵 值 函数 ， 适 当 的 调整 a 的 大 小 可 以 得 到 更 好 的 信号 去 噪 效果。 


3 脉冲 星 信号 去 噪 分 析 


本 实验 中 采用 的 脉冲 信号 来 源 于 Parkes 脉冲 星 实际 观测 数据 ， 该 数据 经 过 PSESTO 去 
干扰 、 消 色散 以 后 得 到 .dat 数据 文件 。 该 数据 文件 以 二 进 制 形式 记 录 ， 采 样 频率 为 1024Hz。 
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图 1 消 色 散 后 得 到 的 信号 图 2 消 色 散 后 得 到 信号 的 功率 谱 


Fig 1 The dispersion 一 limited signal Fig2 The power Spectrum of The dispersion 一 limited Signal 
图 1 中 为 原始 脉冲 信号 消 色 散 以 后 得 到 的 图 形 , 图 2 位 原始 脉冲 信号 消 色 散 后 所 得 信号 

的 功率 谱 ， 本 文 实验 室 对 目标 脉冲 星 信号 文件 按照 Symmlet5 进行 小 波 包 分 解 ， 分 解 层 数 为 
5]. AARP SD WITTER. PRE. POET Bee ae, KH EE BUT sqtwolog 
EHEER WEEET BS EL Ab FE e 
为 了 更 好 的 比较 采用 不 同 的 去 噪 方法 获得 的 信号 去 品 效 果 , 本 文 参照 文献 [13] 采 用 信 品 
IE (SNR), IRERE (PSNR )、 平 滑 度 指标 GD 作为 评价 指标 进行 比较 最 终 产 生 的 去 噪 效 


果 ， 信 噪 比 越 大 、 峰 值 信 噪 比 越 大 、 平 滑 度 指 标 越 小 则 去 噪 效果 越 佳 。 具 体 公式 定义 如 下 : 
(1) 信 噪 比 指 的 是 信号 和 噪声 的 具体 计算 公式 为 ; 

SNR = 10 logo [E2 f(D?) /C1 s? — FO] (3-1) 
(2) 均 方 根 误差 的 具体 计算 公式 为 ; 


RMSE = E m (s(i) — f (D)? (3-2) 


(3) 峰值 信 品 比 指 的 是 信号 的 最 大 功率 和 噪声 功率 之 间 的 比值 ， 用 PSNR 表示 ， 具 体 公式 


PSNR = 10log,o[max(f ())?) / -« (L218)? — f()?)] (3-3) 


(4) 平滑 度 指 的 是 去 噪 后 信号 的 差分 数 的 方差 根 和 信号 之 间 的 差分 数 的 方差 跟 之 间 的 比 
值 ， 常 用 r 表 示 ， 具 体 公式 为 ; 

XP G+ D - FOIE + D-sO)?7] GD 
以 上 计算 公式 中 s( 站 表示 含 噪 信号，F(D) 为 去 噪 后 的 信号 。 
以 下 分 别 为 小 波 变换 去 噪 法 、 小 波 包 变 换 采 用 不 同 阔 值 处 理 法 处 理 信号 后 的 结果 图 、 功 
率 谱 图 以 及 去 噪 后 的 各 评价 指标 汇总 表 。 
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图 3 小 波 软 阔 值 去 噪 后 的 信和 号 
Fig3 the result of de-noised with soft threshold 
= 
图 4 ssc Ti Be (Ss tS 
Fig 4 the result of de-noised with wavelet soft threshold 
1o41 symietsS 
= 


图 5 小 波 软 、 硬 阔 值 折 中 法 去 噪 后 的 信和 号 


Fig 5 the result of de-noised with wavelet soft and hard threshold 
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Fig6 power spectrum of the result of de-noised with wavelet packet threshold method 
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图 7 “LK BUSS UR Wf 
Fig7 the result of de-noised with wavelet packet soft threshold 
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图 8 MEERE A R fi S 
Fig 8 the result of de-noised with wavelet packet hard threshold 
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Fig 9 the result of de-noised with wavelet packet soft and hard threshold 
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Fig10 power spectrum of the result of de-noised by wavelet packet threshold method 
表 3.1 小 波 、 小 波 包 变换 去 噪 效果 表 
Table 3.1 De-noised effect table of wavelet and wavelet packet transform 
index r SNR PSNR RMSE 
Wavelet packet soft 
0.000068 128.751565 50.216641 1.1590332361351423e+38 
threshold method 
Wavelet soft threshold 
0.000068 128.568456 50.205368 1.1590069073995974e+38 
method 
= 4 Wavelet packet hard 
— 0.000068 128.574093 50.216641 1.1590332361351423e+38 
threshold method 
- Wavelet hard threshold 
0.000070 127.619843 48.308142 1.1589152710109917e+38 
method 
um, Wavelet packet soft and 
e 0.000068 128.574093 50.216641 1.1590332361351423e+38 
CN hard thresholding method 
Wavelet soft and hard 
0.000069 128.273024 49.614505 1.1589381807874201e438 
centring method 


BE 3.1 HY AID De CB FE TS E Cp E fi 2: 08 TA A E 2 AE ASR AE A 

= 对 较 好 的 效果 , DS AE NR LAER PFS RR V EE SEAR AE TD 

£ 明显 , 但 是 在 信 品 比 、 峰 值 信 品 比 以 及 平滑 度 方面 都 有 相应 的 提升 。 如 实验 显示 ,采用 小 波 
包 硬 阔 值 法 与 小 波 变 换 硬 阔 值 法 相 比 ， 其 平滑 度 降 低 至 0.000068、 信 噪 比 提高 至 
128.574093dB、 峰 值 信 品 比 提高 至 50.216641dB, ， 小 波 包 软 阔 值 、 软 硬 折 中 浆 值 法 去 噪 在 去 
噪 效 果 的 评价 指标 上 都 有 其 相应 的 提升 。 


4 eZ 吉 与 展望 


本 文 针对 小 波 变 换 在 去 噪 方面 存在 的 不 是 ,提出 了 使 用 小 波 包 闵 值 法 对 脉冲 星 信号 去 噪 
进行 处 理 。 小 波 包 阐 值 法 处 理 首先 需要 对 小 波 包 分 解 后 的 系数 进行 按照 从 小 到 大 的 顺序 排 
列 ,之 后 采用 固定 阔 值 准则 求 出 阔 值 ， 利 用 常见 的 软 阔 值 函数 、 便 闪 值 函数 、 软 便 折 中 浆 值 
函数 对 小 波 包 变换 后 的 小 波 包 系数 进行 处 理 , 利用 处 理 后 的 信号 进行 重 构 ,从 而 得 到 去 噪 后 
的 信和 号。 文章 首先 从 理论 方面 对 小 波 包 几 值 法 去 噪 进行 了 相关 论述 , 之 后 通过 实验 进行 验证 
了 小 波 包 闵 值 法 去 除 脉冲 星 信 号 噪声 的 效果 优 于 小 波 分 解法 。 为 了 使 实验 结果 得 到 量化 标 
准 ， 本 文采 用 了 信和 号 去 噪 效 果 评 价 指标 信 噪 比 、 峰 值 信 噪 比 、 平 滑 度 以 及 均 方 根 差 进行 量化 
比较 不 同方 法 去 噪 的 效果 。 实 验 结果 表明 小 波 包 闪 值 法 去 噪 较 小 波 变换 有 更 好 的 取消 效果 。 
由 于 目前 本 文 只 使 用 了 Symmlet5 小 波 基 和 固定 阐 值 准则 对 选用 的 Parkes 脉冲 星 信号 进行 了 


— 


T 


Xl 


测试 , 是 否 适用 于 其 他 脉冲 星 信号 、 是 否 其 他 小 波 基 或 者 闪 值 准则 也 能 够 获得 同样 的 去 噪 交 
M JE. 
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Research on denoising method of pulsar signal based on Wavelet Packet 
Threshold Method 


Shao Guangsheng', Liu Zhijie’!, Yu Qiuyu’, youshanping! 
1 (Guizhou Key Laboratory of information and computing science, Guiyang 550001) 
2 (School of Cyberspace Security, Sichuan University, Chengdu 610207) 


Abstract: Wavelet transform de-noising only decomposes and processes the low-frequency part 
of the signal, neglects the decomposition and processing of the high-frequency part of the signal, 
which can not effectively extract the whole frequency band information of the signal, thus 
affecting the denoising effect of pulsar signal. In view of this, this paper proposes to use wavelet 
packet threshold method to analyze pulse signal, and selects different threshold functions to 
process the coefficients after wavelet packet decomposition. The experimental results show that 
the denoising effect of pulsar signal is improved by using this method, and the signal-to-noise 
ratio, peak signal noise and smoothness index of pulsar signal are improved correspondingly, 


which provides a new idea for the denoising of pulsar signal. 


Keywords: Wavelet packet decomposition; Threshold method; Pulsar signal; Denoising 
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